Decoherence of two-electron spin states in quantum dots 
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The time evolution of spin states of two electrons interacting with a nuclear spin bath in a 
quantum dot system is studied. The hyperfine interaction between the electrons and the nuclear 
spins is modeled by an isotropic Heisenberg interaction, and the interaction between the electron 
spins by Heisenberg exchange. Depending on the extent of the overlap between the spatial wave 
functions of the electrons, there are two physically different cases, namely the two qubits either 
interact with the same set of nuclear spins or they see different nuclear spin environment. In the 
two cases, the decoherence of the two-qubit state is studied analytically. We have identified a class 
of two-qubit states which have a rich dynamics when the exchange interaction between the qubits 
becomes large in comparison to the hyperfine interaction strengths. The decoherence time scale is 
determined as a function of the bath-spin distribution and the polarizations of the initial two-qubit 
state. States with large decoherence times are identified by performing a minimization over all the 
two-qubit pure states. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Decoherence in quantum systems is a major obsta- 
cle for quantum computation and information process- 
ing, causing errors during the implementation of vari- 
ous gate operations. Decoherence is well-studied in sin- 
gle qubit systems which are realized by pseudo two-level 
systems in Superconducting Quantum Interference De- 
vices (SQUIDS) [1], Nuclear Magnetic Resonance (NMR) 
@, and Quantum Dot [j| systems. Although single 
qubits form the fundamental unit for quantum comput- 
ing, multi-qubit systems are required to perform non- 
trivial gate operations and implementing quantum algo- 
rithms [U. The study of decoherence in multi-qubit sys- 
tems is, therefore, quite important. It is a more involved 
task here, as one has to look at the effects of environment 
on both the loss of coherence and the entanglement of the 
initial state. As entanglement is believed to exponentially 
speed up the implementation of quantum algorithms [5| , 
the time scales for the entanglement loss are as important 
as the decoherence time scales. 

In this work we study the spin dynamics of two elec- 
trons confined in a quantum dot geometry. That is, we 
inject a pair of electrons (in a definite spin state) into a 
confined region of nuclear spin bath. The two electrons 
are coupled to each other by an exchange interaction, and 
to the nuclear spins through a contact hyperfine interac- 
tion. The exchange interaction strength J is a measure 
of the overlap of the spatial wave functions of the two 
electrons (qubits). There are two distinct physical cases 
we have to consider. The spatial wave functions of the 
electrons can have a strong overlap, indicating a large 
value of J. Also both electrons interact with the same 
set of neighboring nuclear spins here, which indicates a 
common nuclear bath for the two-qubit system. On the 
other hand, when the two electrons are physically apart, 
the overlap is quite small indicating a smaller value of 
J. Here, each electron interact with a different set of 



neighboring nuclear spins, and thus implying a different 
spin bath for each electron. Since J can be controlled 
by gate voltages one can tune from one regime to the 
other. Recent experiments [y, [7| have shown a strong 
dependence of decoherence time scales on the exchange 
interaction between the qubits. In a semi-classical treat- 
ment given by Taylor et al [8|, the interaction of the 
qubits with the nuclear spins is replaced by an effective 
magnetic field. The decoherence time scales are found 
by performing an ensemble average of the spin expecta- 
tion values over different configurations of the magnetic 
field. A recent review by Hanson et al [9( summarizes the 
results of experimental procedures employed in studying 
the spins in few-electron quantum dots. Gywat et al |l0| 
have studied the evolution of qubits interacting with an 
off- resonant cavity (instead of a nuclear spin bath). De- 
coherence of two electrons, interacting with different nu- 
clear baths, in the presence of a magnetic field has been 
studied by Zhang et al The decay of the singlet- 

triplet correlator in a double dot system was studied by 
Coish and Loss [Hj]. The evolution of a central two- 
qubit system interacting with a spin bath is studied by 
several authors [H, [13] • A bath- induced entanglement 
between initially unentangled qubits has been studied 
earlier [l5|, [l^|. It was shown by Storcz et al [T3| that 
the errors during gate operations involving coupled solid 
state qubits can be reduced when both the qubits are 
coupled to a common bath at low temperatures. Effects 
of decoherence on the transfer of quantum states along a 
spin chain are studied in [l8| , where the authors have con- 
sidered separate spin bath environments for each qubit. 

In our earlier work [l9|, |2C| we have shown that the hy- 
perfine interaction between a qubit and the nuclear spin 
in quantum dots can be mapped to an effective isotropic 
Heisenberg interaction between the qubit and the neigh- 
boring nuclei. Such an approximation is valid over time 
scales where the contribution of the sub-dominant terms 
in the Hamiltonian are weak. Under the same approxi- 
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mation, a model Hamiltonian in the present case of two 
qubits interacting with the nuclear spins can be written 



"H = (K a Sa ■ h A + K B S B ■ Ie B ) + JSa- S[ 



(1) 



where the qubit spins Sa, $>b couple to the total nuclear 
spins Ig A , Is B , with respective coupling strengths Ka 
and Kb- It should be noted that Ig A obtained by sum- 
ming over all the nuclear spins that interact with the 
qubit A through a contact hyperfine interaction. The 
exchange interaction between the two qubits is given by 
the last term in EqfT] As we argued earlier, depending on 
the exchange interaction strength J, we will consider two 
physically different situations, viz. (i) a common bath for 
the qubits and Ig A = Ig B (when J is large, i.e. the spatial 
wave functions of the qubits strongly overlap) (ii) differ- 
ent baths for the two qubits and Ig A Ig B (when J is 
small, i.e. the two qubits are physically far apart). For a 
GaAs quantum dot, typical values of Ka and Kb [1] are 
O(10~ s eV). Depending on whether the bath is common 
or different, the exchange interaction strength J can take 
values from a few to a few lOneV. The number 

of nuclei A with which the qubits are strongly interact- 
ing, can vary too. When the baths are isolated, N ~ 10 3 
2l|. On the other hand, when the overlap between the 
electrons is very large, they together see a larger nuclear 
bath with A ~ 10 5 . 



A. The initial state 

We take the initial qubit-bath state to be a direct prod- 
uct, p(0) = pab(0) ® Ps(0)- Pab can be represented as 



Pab 



X 
4 



1 1 3 

P a . Sa + -Pb-Sb+ ]T IL mn S%S B , (2) 



where the vector polarizations are given by Pa,b = 
2Tt[pabSa,b] 1 and the cartesian components of the ten- 
sor polarization are II™" 1 = 4Tr[pABS A S B ]. We note 
that for a pure state {p 2 AB = Pab), we have Pa = Pb < 1, 
and Pi + P§ + EmJn™") 2 = 3 - :t should be noted 
that pure states have a non-vanishing tensor polariza- 
tion. If pab is also maximally entangled, then we 
havePA = Pb = 0, implying maximal tensor polarization 
strengths. On the other hand, for direct product pure 
states, we have Yl mn = P^P^. p AB for a maximally- 
entangled state p E AB (Q) =%[\ U) - | |T}][(U I - (IT I], 
and a direct product state Pa B (0) = | T-lXt-l I can be 
written in terms of the various polarizations as 



Pab 



u 
Pab 



- A [I-A{SIS% + S%SI 



I 

4 

\[1+2{Sa-S b )-AS a S b ]. 



(3) 



We now model the initial state of the bath follow- 
ing the approach taken by Rao et al We write 



the state of the bath as an incoherent superposition of 
states labeled by the bath spin Ig , with weights Xi £ , 
PS (0) = Ph (0) ■ I 11 this study all pj £ (0) will be 

taken to be unpolarized (multiple of identity). The 
weights Xi £ are however, free parameters. This is not too- 
restrictive an assumption since we have shown in [l9| that 
higher-order polarizations in each spin sector of the bath 
give small corrections to the decoherence time scales. De- 
pending on whether the bath spins interact ferromagnet- 
ically or antiferromagnetically (j|) can be large (~ A 2 ) 
or small (close to zero) , where A is the number of nuclear 
spins. For a completely unpolarized bath with pg = iprl, 
the bath-spin distribution is Xi £ w ij exp(— i|/2A). For 
this state (If) ~ N. 

The dynamical evolution of the system is governed by 
the equation p(t) = Up(Q)U\ where U is the time evo- 
lution operator. The two-qubit reduced density matrix 
Pab {t) is obtained through a partial trace over the bath 
degrees of freedom, pab {t) can in general be a mixed 
state, which would then imply that Ti p AB (t) < 1. We 
can use the extent of mixing as a measure of decoherence 
which is given by 



D(t) = l-Tvp 2 AB (t) 
= ~[3-(P2(t)- 



mi i 



(ir™) 2 (t))]. (4) 



The above measure has a minimum value of zero for a 
pure state and a maximum value of | for the completely 
mixed state (pab = As a measure of entanglement, 
we shall use the concurrence measure [22] in our analysis. 
For a pure state, the concurrence is simply given by C — 
1 — P\. Even though the expression for the concurrence 
for a general mixed state is not that straightforward, it 
can be expressed in terms of the P\ b an d W mn . For 
example if [pabjS\ + S B ] = 0, then the concurrence is 
given by 



C = —max 
2 



Uw) 



(n 



xy 



x\2 



IP*)2 + (P 1+P J 



(5) 



In the next section, we study the dynamics for non- 
interacting qubits. We determine the evolution of the 
the reduced density matrix pab analytically and thereby 
extract the decoherence measure D(t). In Sec-III we con- 
sider the case of interacting qubits (J ^ 0). The dy- 
namics becomes more complicated, and is sensitive to 
whether the couplings of the qubits with the bath are 
same or different. We study both these cases, and the re- 
sults are again obtained analytically. A short-time anal- 
ysis of decoherence is made, which allows us to identify 
a class of states with large decoherence time scales. We 
evaluate the decoherence measure and the concurrence 
for longer times in Sec-IIIA (Ka = Kb) and in Sec-IIIB 
(Ka^Kb). 
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II. DYNAMICS WITH NON INTERACTING 
QUBITS 

In this section we consider the case of non interacting 
qubits, i.e., the overlap between the spatial wave func- 
tions of the qubits is negligible. Thus, we set J = 0. 
As argued in Sec. I, the nuclear environments for the two 
qubits are different (Ie A ^ Is B ) ■ The Hamiltonian is 

H = K A S A ■ I £a + K B S B ■ h B ■ (6) 

Since the two qubits evolve independent of each other, 
an initial direct product state will evolve into a direct 
product state for later times, i.e. if pab{0) = Pa(0) <8> 
Pb(0), then we have pab (t) — pA{t) ® Pb{P). Hence 
the dynamics is of interest only when the initial state is 
entangled, which we study below. 



Following the method developed in |l9j, we can write 
down the unitary operator, that evolves the qubit-bath 
system, as 

U (t) = (pa + q A S A -h A )® (pb +q B S B - h B ). (7) 

Here the timc-dcpcnclcnt coefficients are given by pa (t) = 
cos A^i + iKa sin h. A t/2kA, q A {t) = 2ii^^sin A^t/A^, 
where 2Aa = Ka{Is a + 1/2), and a similar form for 
PB(t),qB(t) of the qubit B. Using the above time- 
evolution operator, we evaluate the time dependent po- 
larizations of the two-qubit state, and we obtain 

P A {t) = gi(t)P A (0), P B (t) =gi(t)P B (0) 
Tl mn (t) = g 2 (t)U mn (0). (8) 

The time-dependent coefficients are given by 
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(*) = E A feA A feB | PAP B| a {l- f£a(J S + 1) 

Ie A Js B I 
(*) = E ^ A ^>AP B Al- l£B{h ° +l) 
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PAPB 





(9) 



The coefficient g\ (t) is obtained from gi (t) by interchang- 
ing the labels A and B. We note that g\(t) > g%(t) 
i.e., the initial tensor polarization n mn decays faster in 
comparison to the vector polarization. The decoherence 
measure (given in Eq[4]) for the present case has a simple 
form and is given by 

D(t) = ~{3(1 - gl{t)) 2[gl(t) - g 2 2 {t)]P{Q) 2 }. (10) 

where P(0) = P A (0) = Pb(0). Since the initial con- 
currence is given by C(0) = 1 — P(Q) 2 , it is clear that 
decoherence is sensitive to the entanglement in the ini- 
tial state. Decoherence becomes stronger with increas- 
ing initial entanglement. These results are shown in 
Fig.l where we have plotted D(t) for states with differ- 
ent initial entanglement. We can now extract the time 
scales for decoherence and the entanglement loss from 
the short-time behavior of the above time-dependent co- 
efficients. For this, we consider Ka — Kb — K, and 
also = (Is B ) = (Is) i- c -; the individual baths have 

the same bath-spin distribution, and each qubit is inter- 
acting with its bath with the same interaction strength. 
The leading-order time dependence of the coefficients are 
given by 

gi(t) - i~\K 2 (il)t 2 

g 2 (t) « 1-|K 2 (/I)i 2 (11) 



The decoherence measure for short times is then given by 
D(t) ps 1 — exp— £ 2 /t!,, with the decoherence time scale 
given by 

\ = l K 2 (ii)(:z-p(vY) (12) 

As can be seen from the above, any entanglement 
present in the initial two-qubit state would imply P(0) < 
1, which in turn will decrease the decoherence time scale. 
Thus, we conclude that initial separable two-qubit states 
decoher less compared to the entangled states when the 
exchange interaction between the qubits is negligible, and 
the two qubits interact with different sets of nuclear spins. 
It should be noted that the decoherence time scale de- 
pends on the initial state of the bath and the qubits, 
which is a hallmark of a non-Markovian evolution p^ |. 
For initial states with a nonzero entanglement, i.e. the 
concurrence at any later time has a complicated struc- 
ture. 

However, for an initial two-qubit state | jj, +r J,|> 
I \J\ + r 2 , which is an eigenstate of S z , the concurrence 
is straightforward to calculate using Eqs.5, 8 and 9. The 
initial polarizations here are given by P(0) = (1— r 2 )/(l+ 
r 2 ),W z = -1,U XX = W y = 2r/(l + r 2 ). The short- 
time behavior of the concurrence, using the short-time 
expansions of the time-dependent coefficients g\ and <?2, 
is C(t) w C(Q) exp(—t 2 /tq), where the decay time scale 
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FIG. 1: Purity measure Tr[p AB (t)] for a state with an initial 
concurrence C(0), is plotted as a function of time. The initial 
state of the baths are unpolarized consisting of N = 100 spin 
1/2 nuclei. It can be seen from the figure that, the loss of 
purity is faster for entangled states. Also plotted here is the 
concurrence C(t), for a maximally entangled state. The loss of 
entanglement is faster than the loss of purity for a maximally 
entangled state. We have taken Ka = Kb = K. 



is given by 



2P(0) 5 



3^VI) 1 _ P{()) , 



It should be noted that for initial unentangled states, the 
concurrence is zero for all times, as the evolution of the 
two qubits is independent of each other, and the state is 
a direct product for all times. It is clear from the above 
expressions for the time scales that r c < 775 for the en- 
tangled states, and the equality holds for the maximally- 
entangled states (P(0)=0 here). Even though C(t) and 
D(t) have a similar short-time behavior, the concurrence 
falls off rapidly to zero for later times, as shown in Fig.l. 
The entanglement in the state vanishes when g2[t) < 1/3. 



III. DYNAMICS WITH INTERACTING QUBITS 

In this section we consider the case of interacting 
qubits. This corresponds to a considerable overlap of 
the spatial wave functions of the two qubits, implying 
a nonzero exchange interaction, and consequently the 
qubits see a common nuclear bath (I £a = Is B = Is)- 
The exchange interaction is larger than the Overhauser 

held i.e., J > (K A + K B )\J (If). Thus we write the 
Hamiltonian as, 



The time-evolution operator in this case can be written 
as, 



U = 



C2 



ai(t)+a 2 {t){S A -SB)-h\ (1- 
03 (t) + ai(t)SAB ■ h + a 5 (t)(SAB ■ h] 2 

a 6 (t)(S A - S B ) ■ I £ + a 7 (t)(S A x S B ) ■ I £ 



C2 

2 

(15) 



The expressions for the coefficients ai(t) are given in Ap- 
pendix A (Eq |A8[) . Note that the coefficients a2(t),a,Q(t) 
and aj (t) survive only when there is an asymmetry in 
the qubit bath coupling strengths (K A ^ Kb)- The op- 
erators (S A — Sb) ■ Is and (<SU x <9b) * Is appearing in 
Eq|15l cause transitions between the singlet and triplet 
subspace of the qubits. 

The time dependent polarizations are given by 



Pa® 



h(t)P A (0) + f 2 (t)P B (0) 
3 3 

+/ 3 (t)$>™ E e n«n«(o), 

n—1 — 1 

P B (t) = U(t)P B (o) + h(t)P A (o) 
3 3 

~U(t)J2^ E CnijWiO), 

(13) u mn (t) - / 7 (t)n™(o) + / 8 (t)n™(o), 

+/ 9 (i)5 mn Trn + f w {t)e mnk {P%{Q) - P|(0)). 

(16) 



where e„ stands for a cartesian unit vector. The ex- 
pressions for the coefficients fi(t) are given in Appendix 
IB"1 We extract the decoherence time scale td from the 
leading-order time dependence of the coefficients fi (t) for 
short times, and we have 

h{t)Kl-{{P £ )Kl/Z + J 2 /i)t\ 

f2(t)J 5 (t)^J 2 t 2 /A, 

h(t),Mt),f 10 (t)^Jt/2, 



kit) 

hit) 
hit) 
hit) 



\-{{i 2 £ )Klii + j 2 /A)t\ 

l-{{ll){K\+Kl)/Z + J 2 /A}t\ 

-((I 2 )K A K B /3-J 2 /4)t 2 , 

{I 2 )K A K B /3t 2 . 



(17) 



Using the above short-time behaviour of the coefficients, 
— ( — — ) 2 

we have D(t) ~ 1 — e T o> f and the decoherence time 
scale td is given by 



- - Up) 



[(K 2 A + K 2 B )(3-P 2 (0j) 



H = (K A S A + K B S B ) ■ h + JS A ■ S B - 



(14) 



-^^ B (Trn 2 (0) - (TrII(0)) 2 )] . (18) 
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Let us recall that P(0) = P A (0) = P B (0) for pure two- 
qubit states. If the initial state is separable, then we 
further have, TrII 2 (0) = (TrII(0)) 2 . The decoherence 
time scale for these states is given by 

(i) =\(i 2 eKK A +Kl)^± + ±, (19) 

\ T D/sep 6 T A T B 

where ta and tb represent the decoherence time scales 
for the individual qubits. Clearly the individual qubits 
decoher slower in comparison to the two-qubit separable 
state. For maximally-entangled initial states (P(0) = 0), 
the decoherence time scale depends on 1Z = Trll 2 (0) — 
(TrII(0)) 2 , which further depends on the details of the 
state. For example, 1Z = 2 for the Bell states in the 
triplet sector, and 1Z = —6 for the singlet state. For the 
two cases, the decoherence time scales are respectively 
given by, 

^) S = -2 {P£)(KA - KB)2 > 

(4) = \^)( k a + k% + \k a k b ). 

(20) 

The time scales are sensitive to the relative sign between 
K A and Kb- If K A and Kb both have a positive sign (an- 
tiferromagnetic interaction) or a negative sign (ferromag- 
netic interaction), then the triplet sector suffers a faster 
decoherence, while the singlet state decohers faster if they 
have opposite signs. It can also be seen from Eq [TO1 and 
Eq|20l that the separable states decoher slower than the 
entangled states if K A and Kb have opposite signs. For 
partially-entangled states, the calculation becomes more 
involved, and a global minimization over the initial two- 
qubit pure states for finding the least decohering states 
will be done in Sec-IIID. For the maximally- entangled 
Bell states (C(0) = 1), the concurrence for later times 
is given by C{t) = max{ \W x {t)\ - |1 + IF z (i)|/2, 0}. It 
has a small-time behavior C(t) w e~^^ Tc \ with a de- 
cay time scale equal to the decoherence time scale given 
above for the singlet and triplet state. Thus, the loss of 
entanglement and decoherence occur over the same time 
scale. This feature for the maximally- entangled states, 
that tc = tb, we have seen in the case of the qubits 
interacting with different nuclear baths also, except that 
there all the bell states have the same decoherence time 
scale. Below we examine the dynamics for larger times, 
t > tb in detail, with symmetric and asymmetric cou- 
plings separately. 



A. Identical coupling strengths K A = Kb ~ K 

When K A = Kb, the dynamics becomes simple, as the 
triplet and the singlet components of a two-qubit state 
evolve independently. The singlet sector is unaffected, 



and a nontrivial evolution occurs only in the triplet sub- 
space. The time-evolution operator given in Eq|15l sim- 
plifies now to 

u = (i_5k) + [ a3 ( <) + a4 (i)^ B . / ; 

+a 5 (t)(S AB -f £ ) 2 }^. (21) 

We now consider the evolution of the initial two-qubit 
states either maximally-entangled or unentangled sepa- 
rately. 

1. Unentangled states 

For the unentangled state, let pab(Q) = I T4)(t-l I- The 
initial polarizations of this state are P\ — —P B — 1 an d 
IF 2 = — 1. With time, the system picks up other compo- 
nents of tensor polarization. The nonzero polarizations 
and the concurrence at any later time t are given by 

Pl(t) = -PUt) = [fi(t)-f 2 (t)]Pl(o) 
n zz («) - [Mt) + / 8 (t) + /9(t)]n"(o) 
n xx (t) = u yy (t) = f 9 (t)n zz (o) 

W y {t) = _n»*(t) = / 10 (t)[P*(0) - P B (0)] 
C(t) = ^max^ fi(t) + fUt) 

-|l-/ 7 (t)-/ 8 (t)-/ 9 (t)|, |, (22) 

where the coefficients fi(t) are given in Appendix B. In 
Fig. 2 we have plotted the polarizations P A {t), H xx (t) 
and H zz (t) as a function of time. We have also plot- 
ted the concurrence in the same figure. One can see that 
P| shows an under-damped behavior with a saturation 
value oscillating between The oscillation frequency 
is 2ir/J. The behavior of the off-diagonal component of 
the tensor polarization H xy is similar to that of P\, but 
phase shifted by 7r/2J. The diagonal components of the 
tensor polarization contribute to the concurrence, but 
have no dependence on J. The initial rapid oscillations 
in the concurrence are due to Tl xy . Hence, even when 
J = 0, the two qubits get entangled. This is an example 
of bath-induced entanglement. 

Since the initial two qubit-state is separable, the deco- 
herence time scale can be calculated using EqQjJ] Taking 
the bath-spin distribution to be A/ £ w 7| exp(— 2/J/7V), 

we get tb = j^y^ - As the individual qubits are also 
initially pure their corresponding decoherence time scale 
is given by t A: b = jtt/^ - As argued above, the two- 
qubit state decohers faster in comparison to the individ- 
ual qubits. Even though the two-qubit decoherence is 
unaffected by the exchange interaction, the concurrence 
and the polarizations of individual qubits are sensitive 
to J. A comparison between the single-qubit and the 
two-qubit decoherence measure is made in Fig. 3 
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FIG. 2: The vector and tensor polarizations are plotted as 
a function of time for an unentangled state interacting with 
an initially unpolarized bath consisting of N = 100 spin-1/2 
nuclei. The time variation of concurrence C(t) is also plotted. 



FIG. 4: The non-vanishing components of tensor polariza- 
tions, the concurrence C(t) and decoherence measure D(t) 
are plotted as a function of time. The initial state of the two 
qubits is maximally entangled, and the bath is unpolarized. 
The bath consists of N = 100 spin-1/2 nuclei. 
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FIG. 3: A comparison between the decoherence measure of 
the two qubit state D(t), and that of the individual qubit 
states DA(t) = Ds(i) = 1 — Tr[p^(t)], is shown as a function 
of time. The initial state of the two qubits is unentangled, 
implying that the individual qubits are in a pure state ini- 
tially. The initial state of the bath is unpolarized consisting 
of N = 100 spin 1/2 nuclei. 



2. Maximally- entangled states 

We shall now consider a maximally-entangled Bell 
state for the two-qubit state, viz. the triplet state with 
S z = 0. The initial polarizations of this state are 
Y[xx _ Y[yy = \ anc j rj zz = — 1. The non- vanishing po- 
larizations and the concurrence at any later time t are 



given by 

n-( t) = ipv (t ) = [ fr (t) + /sW + / 9 (i)]ir-(o) 
n"(t) = [/ 7 (t) + / 8 W-/9W]n Z2 (o) 

C(t) = ^maa:|2|n !BS (t)|-|l-n**(t)|,o| (23) 

where fi(t) are given in Appendix B. In Fig. 4 we have 
plotted the components of tensor polarization, the con- 
currence and D(t) as functions. Setting Ka = Kb = K 
in Eq[5Dl and taking the bath-spin distribution to be 
Aj £ w l\ exp(— 2J|/iV), we have the decoherence time 
scale tjj = j^m- For longer times C(t) is zero, and at 
the instant C(t) becomes zero, the components of the 
tensor polarizations take values fl" = Il aa = 1/2 and 
U zz = 0. 

In the case of non interacting qubits we have seen that 
all the components of the tensor polarizations had the 
same decay time scale. In the present case the trans- 
verse components of tensor polarization (U XX ,H VV ) have 
different decay rates in comparison to the longitudinal 
component Tl zz . This can be seen from the small-time 
expansion of these polarizations which are given by 

n H (() = n yy (t) = [i - 2K 2 {i 2 £ )t 2 ]n xx {Q) 

Il zz (t) = [l-4K 2 (ll)t 2 ]Tl zz (0). (24) 

Denoting the longitudinal and transverse decay rates by 
Ti and T2, we have T\ = 2t 2 . In contrast for the case of 
the noninteracting qubits, we have t\ = t-i. 
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B. Asymmetric couplings Ka 7^ Kb 

In this section we consider the interaction (Ea fTij) in 
its full generality. Here, the qubit-bath interaction term 
does not commute with the exchange term. This causes 
transitions between the singlet and triplet subspaces, due 
to which the singlet state also decohers. In Sec. Ill, we 
have seen that short-time dynamics is independent of J, 
for example, the decoherence time scale does not depend 
on J. However, the effects of a nonzero J become domi- 
nant for times t > l/J, which can only be seen from the 
full solution of the dynamics. For J large, the singlet and 
the triplet sectors of the qubits are effectively decoupled 
from each other. Since the hyperfine coupling strengths 
are small compared to J, the transition matrix elements 
between the singlet and the triplet sectors are negligible. 
Thus, the singlet state and the states with a large singlet 
fraction remain invariant through time evolution. On the 
other hand, the triplet states suffer decoherence, due to 
transitions within the triplet sector. Hence, the situa- 
tions where a triplet is a ground state or a metastable 
state are not favorable for quantum computations. 

We consider a one-parameter class of initial two-qubit 
states of the form 

\1>ab) = 2 * 2) [(1 + r)\So) + (1 - r)|r >], (25) 

in terms of the singlet and triplet (with S z — 0) Bell 
states, i.e. \S ) = U> - | IT)] and \T ) = -^[| U) + 
I ||)]. The evolution of the system is best studied in the 
Bell-basis, which apart from |So) and \Tq) also comprises 
of, |7i) = -y ||) + I ||)], |T 2 ) = -^[| II) - I II)]. The 

two-qubit reduced density matrix state at any later time 
is given by 

p AB (t) = ci(t)|5o)(So| + C3(t)|2b)<To| + ps(*)To)<So| 
+ CittHlTiXTxI + lTaXTaD + CBWlTiXTal 
+ H.C. (26) 

The time-dependent coefficients are given by, 

ci(<) = ^TT^5: A/£[(i+r)2|ai|2+(i - r)2 

h{h + l)(|a 6 | 2 + M 2 + 2Im a*a 7 ], 

C2(<) = ^TT^y^ A/£[(1+r)2(1 ^ |ai|2)/3 

+ (l-r) 2 {|a 3 | 2 + i/ £ (/ £ + l) 
((8/| + 8/ £ -l)^ + 2Re a 3 a* 5 )}}, 



1 _ r 2 2 
C3(*) = 2(1 + r 2 ) ^ Xl£ ^^ a3 + 3 /f( - /g + 1 ' )a5 - > 

+-Is{h + l){a 2 {al+ia* 7 )], 

a(t) = ~[i -<*(*) -<*(*)], 

<%(*) = 6(1 + r 2) Z X ieh{h + l)Re[a* 2 (a 5 + 2a 4 )], 

(27) 

where the coefficients a% carry time dependence, and are 
given in the Appendix A. If the initial state is one of 
the Bell states, pab (t) will be diagonal in the Bell-basis 
for all times, because the coefficients c^(t) = c 4 (t) = 0. 
For a general linear superposition state, pab it) picks up 
off-diagonal terms as well. 

To show the effects of J on the states close to the triplet 
and singlet states, we consider two states corresponding 
to r = ±0.5 in Eq.25. The state corresponding to r = 
—0.5, is close to the triplet state, whereas the other state 
with r = 0.5 is close to the singlet. As argued earlier, we 
should expect a large dependence of D(t) on J for the 
state close to the singlet. This, in fact, can be seen from 
Fig. 5, where we have plotted D(t) for these states for two 
values of J. The decoherence measure D(t) for the case of 
r = 0.5, is suppressed by a strong exchange interaction, 
whereas for the former case (r = —0.5), there is only a 
small effect of J on D(t). We note that for r — * 1 (singlet 
state) D(t) —* for large values of J. In the other limit 
when r — > — 1 (triplet state), decoherence persists and 
remains unaffected by the value of J. 

As we have seen earlier in Sec-III, that up to 0(t 2 ) 
there is no dependence of D(t) on J. But when J is 
large the contributions from higher-order terms also be- 
come significant over time scales r < 1/K. The initial 
Gaussian decay is seen only for times t < l/J, and the 
behavior can change drastically for later times. It can be 
seen from Fig. 5, that the behavior of D(t) for the state 
corresponding to r = 0.5 is similar (Gaussian) initially 
(t < l/J) and then changes into a slowly varying oscilla- 
tory function. 

To understand the behavior of D(t) for large J, we 
consider the singlet state, whose decoherence is highly 
sensitive to J. As can be seen from Eg |2"61, the time 
evolved state has a simple structure, we have pab (t) — 
cx (t)\S Q ) (So I + (1 - Cl (t))[\To) (To \ + \T 1 ){Tx\ + \T 2 ) (T a |]/3, 
where C\(i) is given by 

d(t) = ^2 A/ £ [cos 2 A_t +p 2 sin 2 A_i]. (28) 
le 

The coefficients A_ and p are given in Appendix A. For 
J ^> \/N(Ka + Kb), we can approximate the coefficients 
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FIG. 5: A plot of D(t) for the partially entangled states corre- 
sponding to r = ±0.5 (see Eq |25|l is shown for different values 
of the exchange interaction strength J. The initial state of 
the bath is unpolarized consisting of TV = 100 spin 1 /2 parti- 
cles. Here K — (Ka + Kb)/2 is the average of the qubit bath 
interaction strengths. 



A_ and p as 



A_ 



(K A -K B f 2 

P ~ 1-2 2J2 I £ 

For a bath-spin distribution given by Aj £ 
J|exp(-2/|/7V), we find 



(29) 



Cl (t) = l - 3/3 2 



1 



cos{Jt + § tan _1 /3i) 
(l+/3 2 i 2 )t 



(30) 



where /3 = (Ka — Kb)V N/2J. Here the decoherence 
measure is given by D(t) = 2/3(1 — c\(t)). For large J, 
since Ci(t) ~ 1, -D(i) is quite small, implying that the 
singlet state hardly decohers. As can be seen from the 
above, there is a slow- varying oscillatory behavior for the 
decoherence measure. A similar power-law behavior for 
the singlet-triplet correlator for large J was shown earlier 
by Coish and Loss [l2| . This also explains the oscillatory 
behavior of D(i) for large J for the state corresponding 
to r = 0.5 (which is close to the singlet) seen in Fig 5. 

The Bell states, which are widely used in quantum 
information protocols, rather have a special evolution. 
As noted earlier these states remain diagonal in the Bell 
basis for later times. For a Bell state we have Pa — 
Pb = 0, and only the diagonal components of n to be 
non vanishing. The state at any later time has the form 



1 



PAB(t) = j[i + ^n mm (t)s^s^]. 



(31) 



It follows from EqOH that W nm (t) = [f T (t) + 
f 8 {t)]U mm (0) + f 9 (t)TvU. This is true only when the 



bath is unpolarized. The state can pick up other polar- 
izations if there is a nonzero magnetic field, or when the 
bath is polarized. We do not consider such cases in this 
work. 

Before we end this section we consider a class of 
maximally- entangled mixed states, viz. the Werner 
states. A Werner state is an admixture of the singlet 
and the identity state viz., pab = p\Sq) (Sq\ + —jrl- 
where |5o) is the two-qubit singlet state. This state re- 
mains invariant under time evolution for all values of p 
(0 < p < 1). The concurrence for the above given state 
C = max[(3p — l)/2, 0]. Werner states are entangled for 
p > 1/3 and are separable for p < 1/3. When Ka = Kb 
the Werner state remains unchanged as the singlet has no 
evolution. However, when Ka ^ Kb, the state evolves 
as the singlet decohers. The bath effects on this state 



can drastically get reduced when J ^> (Ka + Kb) 



(it) 



C. Least-decohered two-qubit pure states 

We now address the problem of identifying the least 
decohering two-qubit states through the time evolution. 
The initial state of the system is a direct product of the 
two-qubit pure state (which may or may not be a direct 
product) and the bath state (which is a direct sum over 
bath states from various spin sectors with a bath-spin 
distribution). We make no assumption on the relative 
magnitudes or signs of Ka, Kb and J. The general ex- 
pression for the decoherence time has been obtained ear- 
lier from the short-time expansions of the time-dependent 
coefficients (see Eqn lT5|) . The decoherence time scale for 
a two-qubit state is given by 



1 



\{P £ ){({KaSa + K b Sb) 



(K a S a 



K b Sb) 2 }, 
(32) 



where the expectation values are evaluated in the initial 
direct product state of the qubits and the bath. It should 
be noted that the decoherence time scale is independent 
of J, as the leading-order contribution for small t is in- 
dependent of J for all the time-dependent coefficients. 

We wish to minimize the above expression over all the 
pure states. To that end, we write the most general pure 
two-qubit state as 



1 



(33) 



where h\, fi2 are two arbitrary directions in space. Here, 
7 is any complex number which determines the entan- 
glement. The ket | f^ t ) is an eigenstate of the operator 
Si.hi with an eigenvalue +1/2. Since the Hamiltonian 
is isotropic, and the bath is unpolarized, we can choose 
fii = z, without any loss of generality. Let us denote 
the polar angles of the other arbitrary unit vector fi2 by 
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(6, 4>). The decoherence time scale for the above state is 
given by 



1 



I 2 £ ){K\ + Kl) 



2|7| : 



(l + M 



2\2 



(1 — 8 cost 



-2<5cos 2 (0/2) 



i + M 2 



(34) 



where we have introduced a hyperfine-inhomogeneity in- 
teraction parameter 8 — M^-. The decoherence time 

scale is independent of (f>, and it it is easy to check that 
the optimal value of the decoherence time is for 9 = 0, 
i.e. n 2 = z in the above state. It is also straightforward 
to check that 7 is real for the least-decohered state. This 
would imply that it suffices to optimize the decoherence 
time for a state of the form \tp >= | fj, —7 |t)/(l + 7 2 )> 
with 7 real. From the examples considered in the earlier 
sections, one would have expected the above form for the 
least decohering state. We can choose the value of 7 by 
optimization of the decoherence time scale as a function 
of the coupling constants. In Sec-IIIC we saw that states 
with 7 > (see EqfJH]) have long decoherence times when 
K a ,Kb are of the same sign i.e., 8 > 0. Maximizing the 
decoherence time scale, we find the optimal value of 7 to 
be. 



(1-8)- VT=28 



Topt 



7o P t = l; 2 - § - 1 - 



-i<*4 



(35) 



For 8 > 1/2, the least-decohered state is the singlet 
state (7 = 1). For 8 = 0, which corresponds to either 
Ka or Kb is zero, i.e. one of the qubits is decoupled 
from bath, the direct product states decoher slower than 
the entangled states. In Fig. 6, we have plotted the de- 
coherence time for the maximally-entangled states (both 
for the singlet with 7 = 1 and the triplet with 7 = — 1) 
and the separable state (7 = 0) as 8 is varied. Also, 
we have plotted the best possible decoherence time (note 
that the state itself is determined by jopt)- The inset in 
Fig. 6 shows the variation of 7 op t as a function of 8. It can 
be seen from the figure that separable states have a larger 
decoherence time compared to the maximally-entangled 
states over a large range, i.e. 8 < 1/3, however, the state 
with the largest decoherence time has an intermediate 
value of entanglement. 

The above analysis can be generalized for a more gen- 
eral Hamiltonian 



n = T\k\s a 



K B S B ] ■ I< + JS A ■ S £ 



(36) 



where the interaction strength between the qubit A(B) 
and i'th nuclear-bath spin is denoted by Ka{B) 1 . The 
short-time analysis can easily be performed in this case 
aslo. The decoherence time scale for an initial two-qubit 
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FIG. 6: Variation of 1/r 2 with 8, for the separable state (7 
; . the singlet (7 = 1) and the triplet (7 = —1) is shown 
above. Also shown is the dependence of 1/r 2 on 8, for the least 
decohering state (Optimal). The values of 7 corresponding to 
the optimal curve are shown in the inset. The square of the 
inverse decoherence time (1/t 2 ) is plotted in units K 2 (Ig}/3. 



state |V>>=|_U 7 7 IT>/(1- 
nuclear spins, is given by 



- I7I 2 ), interacting with N 



£<3 



N 



2|7l 



(1 + I7I 2 ) 2 



(1-A) 



-2A 



Re( 7 ) 



where A 



K> 



(37) 



the 



hyperfine-inhomogeneity parameter. We have introduced 
two scale factors here 771 and 772 (both of which are of or- 
der unity), which depend on the individual interactions 
strengths and also the bath-spin distribution. In general 
it is difficult to calculate the scale factors. For an initial 
state of the bath of N nuclear spins given by ps = ipri, 
it is easy to calculate the scale factors, and we have 
Vi = T /2 = 1- The value of 7 opt for the least-decohered 
state is given by Eq.35, except that the inhomogeneity 
parameter 8 is replaced by A. Taking the electron spa- 
cial wave functions to be Gaussian functions, i.e. the 
groundstate wave functions of a 2-dimensional confining 
harmonic oscillator potential appropriate to the quantum 
dots, we find that the hyperfine-inhomogeneity parame- 
ter to be A w 0.6, when the distance between the qubits 
is of the order of the harmonic confinement length. Thus, 
the case of a common nuclear bath considered in Sec-III 
holds for A > 0.6, as here J ~ 0.1 meV is larger than the 
hyperfine coupling strengths. For A < 0.6, the common 
set of nuclear spins that interact with the two qubits be- 
comes smaller, and the qubits become isolated as A — »• 0. 
In this case, discussed in Sec-II, the exchange interaction 
is negligible, and the two qubits evolve independent of 
each other, interacting with two different sets of nuclear 
spins. 
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IV. CONCLUSION 

The dynamics of two spin- 1/2 particles interacting 
with a nuclear spin bath in a quantum dot is studied. 
The hyperfine interaction between the electron spin and 
the nuclear spins is modeled by an isotropic Heisenberg 
interaction. The time-dependent polarizations that char- 
acterize the two-electron spin (two-qubit) state are calcu- 
lated as functions of the initial two-qubit state, and the 
nuclear bath-spin distribution. When there is negligi- 
ble exchange interaction between the qubits themselves, 
corresponding to the situation where the two qubits are 
physically apart, it is argued that each qubit sees a dif- 
ferent set of nuclear spins as its environment. Here, the 
separable states (i.e. unentangled two-qubit states) are 
found to have have larger decoherence time scales in com- 
parison to entangled states. On the other hand, when the 
exchange interaction strength is appreciable, we argue 
that the two qubits interact with a common nuclear spin 
bath. In this case, we have considered two cases: the 
two qubits have same hyperfine coupling with the nu- 
clear spins (symmetric coupling) and when the couplings 
are different (asymmetric couplings). For maximally- 
entangled states, the time scales corresponding to the loss 
of entanglement and decoherence are the same, irrespec- 
tive of whether the qubits have a common bath or not. 
The effect of a nonzero exchange interaction strength J 
comes is significant in the case of either there is an asym- 
metric coupling between the qubits and the bath or the 
nuclear environment of the qubits are different. We have 
studied the effect of J only in a limiting case, when the 

J S> (Ka + Kb)\J (/J). For intermediate overlap of the 
electron wave functions, where there are some nuclei in 
common and some different, analytical solution to dy- 
namics is a much more difficult task. In such situations 
we have numerically seen that the dominant effect of J 
arises when the initial state is either the singlet or a state 
having large overlap with the singlet. The least decoher- 
ing two-qubit states are found by minimizing the deco- 
herence time scale over all initial two-qubit pure states, 
for a given hyperfine-inhomogeneity interaction. 



APPENDIX A: THE UNITARY OPERATOR 

For the case of qubits interacting with a common nu- 
clear bath, the Hamiltonian can be written as 



n = {K A s A + k b s b ) ■ i £ + j 



Sai 



(Al) 



where Sab is the total spin of the qubits, which can take 
two values and 1. The total spin of the qubits and 
the bath can take three values F = Ig — 1 , Ig , Ig + 1 . 
Denoting the projection operator by Pj,s A b f° r a sector 
with a given total spin F and the qubit spin Sab, we 



have 



-[{Sab-Is) 2 



(7 £ + l)(2/ £ + l) 1 
+{I £ +2)(Sab-Te)+Is 

+(1-I £ )(Sab ■Te) + h] 



C2 



s 2 



AB 



Iefi 



(1 



h{h + 1) 
D AB \ 



[(Sab ■ hf + Sab ■ h\ 



'AB 

2 ' 
(A2) 



There are two sectors with F = Ig, which can mix if 
there is an asymmetry in the couplings (i.e. Ka 7^ Kb)- 
Since F z is conserved, only states with the same F z can 
mix. In the subspace spanned by the two states with 
F = Ig, the Hamiltonian matrix has the simple form 



H 



:k a -k e 



j_Ka± Es ^ Is{l£ + 
VWe + 1)^^ 



(A3) 



Since the matrix elements do not depend on the value of 
F z , there are only two distinct eigenvalues corresponding 
to the two F = Ig sectors. The eigenvalues are given by 

C± - \{J-{Ka + K b )/2) 

±\y/(J ~ [Ka + K B )/2f + Ig(Ig + 1)(K A - K B ) 2 . 

(A4) 



The eigenvalues corresponding to the sectors with F = 
Ig + 1 , Ig — 1 are given respectively as 



A : = J + I £ {Ka + K b )/2, 

\ 2 = J-(Ig + l)(K A + K B )/2. 



(A5) 



Now, the time-evolution operator in the full space can be 
written as 



U 



e' Alt 4+i,i +c a2t P 7£ -i4 +e 4A +*{cosA_£-isinA_t 
[p(l + 2S AB ■ Ig) + iq ^== ]}(P fe ,i + Pi £ ,i), 



(A6) 



where A± = ((+ ± C-)A p = j^, and q = \Jl-p 2 
Thus, the unitary evolution operator can be written as 



U = 



ai(t)+a 2 (t)(S A -S B )-Ig (1-%) 



+ 



a 3 (t) + a 4 (t)S AB ■ Ig + a 5 (t)(S AB ■ h f 



+ a 6 (t)(S A - S B ) ■ Ig + a 7 (t)(S A x S B ) ■ Ig 



S 2 



AB 



2 

(A7) 
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The time-dependent coefficients are given by 

ai{t) = e lA +*[cos(A_i) - ipsin(A_t)], 

zge iA + 4 . . 
= sin(A_f), 

y/We + 1) 



oa(t) = 



03(f) = aiW + ^^e" 1 '"^'), 



04(f) = 

05 (t) = 

a 6 (f) = 

a r (t) = 



-(j|-l)e^*-o 1 (t)(2J £ + l)] ) 

I s (I e + l)(2I s + l) 1 

+ (I £ + l)e^ t ~a 1 (t)(2I £ + l)] , 

iqe lA+t . , k s 
; sin(A_f), 

5- sin(A_£). 

(J £ (J £ + l))f 



(A8) 



For the simpler case of AT^ = ATe, we have A + = A_, 
p = 1, q = implying that a\(t) — 1 and 02(f) = aeit) = 
07(f) = 0. 



in general a complicated structure. However, they sim- 
plify substantially for the case of Ka = Kb, which we 
give below in terms of the time-dependent coefficients 
o»(f), 



APPENDIX B: THE TIME-DEPENDENT 
COEFFICIENTS f t (t) 

The time-dependent coefficients fi(t) (that determine 
the polarizations for any time t) occurring in Eq.16 have 



/o(f) = ^A Z£ {o 3 + 2o 5 
Ie 

hit) = /4 = jE A ^ 



h{h + 1) 



}• 



03 



[|o 5 | 



Is 



+4Re(a 3 og) /£(j£ 3 +1) ]} + \ Re/ (f), 

/ 5 (f) = / 1 (f)-Re/ (f), 
/ 6 (f) = 2/ 10 (f)=Im/o(f), 



/ 2 (f) = 
/ 3 (f) = 

/7(f) = -7 + 7 2^ A/ ^l a3 l 



[|o 5 | 



a 2 (f) = ae (t) = / 8 (f) = /7(f) - Re/ (f), 

hit) 



-4Re(a 3a g)] /g(/£ 3 +1) } + ^Re/o(f), 



g[l-(M*) + / 8 (*))]. 



(Bl) 
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